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Abstract 

Most existing theoretical studies of momentum transport focus on calculating the 
Reynolds stress based on quasilinear theory, without considering the nonlinear momen¬ 
tum flux-(fyfm||). However, a recent experiment on TORPEX found that the nonlinear 
toroidal momentum flux induced by blobs makes a significant contribution as compared to 
the Reynolds stress [Labit et ah, Pliys. Plasmas 18 , 032308 (2011)]. In this work, the non¬ 
linear parallel momentum flux in strong turbulence is calculated by using three dimensional 
Hasegawa-Mima equation. It is shown that nonlinear diffusivity is smaller than quasilinear 
diffusivity from Reynolds stress. However, the leading order nonlinear residual stress can be 
comparable to the quasilinear residual stress, and so could be important to intrinsic rotation 
in tokamak edge plasmas. A key difference from the quasilinear residual stress is that parallel 
fluctuation spectrum asymmetry is not required for nonlinear residual stress. 


I. INTRODUCTION 

Tokamak plasma rotation and toroidal angular momentum transport have been subjects 
of intensive study due to their important role in reducing turbulent transport as well as in 
stabilizing magnetohydrodynamic (MHD) instability, such as resistive wall modes^. Owing 
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to a wide range of beneficial effects on stability, confinement, and performance of tokamak 
plasmas, much effort has been devoted to understanding the mechanisms underlying the 
change in rotation and how to control it. On one hand, there are a number of known causes 
for the plasma rotation to slow down, such as nonaxisymmetric error fields,- 1 ^ loss of mo¬ 
mentum input as a consequence of Alfven eigenmodes,- and edge localized modes (ELMs).- 
On the other hand, toroidal rotation is driven externally by neutral beam injection (NBI) 
on most present day devices. However, beam injection may be of limited utility in providing 
enough external torque in future reactor such as International Thermonuclear Experimental 
Reactor. One alternative is to take advantage of intrinsic rotation (spontaneous, or self¬ 
generated, in the absence of an external momentum input) which has been widely observed 
under a variety of operating conditions.- Consequently, understanding plasma rotation and 
momentum transport under low external momentum input condition is of major interest. 
The total flux of parallel momentum driven by electrostatic turbulence has the forint 

n r ,|| = {n)(v r u\\) + (U\\){v r h) + {v r hu\\). (1) 

Here, on the right hand side (RHS) of Eq. (JT]) , the first term is the parallel Reynolds stress, 
the second term is convection, due to particle flux, and the third term is the nonlinear flux. 
The Reynolds stress can be further decomposed as 

= -x v ^ + m> + n ?r> ( 2 ) 

where, on the RHS, they are diffusion, pinch term and residual stress, respectively. The 

residual stress is thought to be the origin of the intrinsic rotation which has been intensively 

investigated by using quasilinear theory.- In addition, turbulent acceleration is proposed as 

another possible mechanism for driving intrinsic rotation. 1011 Turbulent acceleration acts as 

a local source or sink, which has different physics from the residual stress. The third term 

on the RHS of Eq. (jT]) , (v r nu\\), represents the nonlinear (as opposed to quasilinear) flux, 

driven by processes such as mode-mode coupling and turbulence spreading.- Most existing 

theoretical works on parallel momentum transport neglect the nonlinear flux, which is less 

understood.- 9 However, the nonlinear flux may also influence the rotation profile, especially 

fl 

at the boundary where the relative fluctuation amplitude can be strong, i.e., since- > 1, 

n 0 

the nonlinear flux cannot be dismissed as small. In this sense, momentum transport theory 
is still not well developed. 
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Recent experimental results on TORPEX showed that blob induced fluctuations are so 
strong that the toroidal flow is transiently reversed, and the associated nonlinear toroidal 
momentum flux can be dominant for some time.— This result suggests that the nonlinear 
flux is no longer negligible. Similar blobs in L-mode^r— and edge localized mode (ELM) 
filaments in H-mode^S - — are observed in tokamaks. In general, the nonlinear flux is of po¬ 
tential relevance in the strongly turbulent edge. Therefore, to fully comprehend momentum 
transport, a theoretical study of the nonlinear momentum flux in strong turbulence seems 
necessary. 

In present work, we calculate the nonlinear parallel momentum flux by using the three 
dimensional Hasegawa-Mima equation— containing the compression of ion parallel velocity 
and the ion parallel momentum equation. For comparison, the parallel Reynolds stress is 
also calculated. We find that the nonlinear diffusivity is small compared to the quasilinear 
diffusivity from the Reynolds stress. However, the dominant nonlinear residual stress can 
be comparable to the quasilnear residual stress with opposite sign, if increasing fluctuation 
intensity profile is used for the residual stress.— This indicates that strong momentum trans¬ 
port induced by blob ejection at edge is important to intrinsic rotation. We also find that 
parallel fluctuation spectrum asymmetry is not necessary for nonlinear residual stress, in 
contrast to the case of quasilinear residual stress. 

The remainder of this paper is organized as follows. Sec. II presents the minimal model 
adopted in this work. The nonlinear momentum flux and its comparison to the Reynolds 
stress are presented in Sec. III. Finally, we summarize our work and discuss the implications 
for momentum transport and rotation response to blob ejection in Sec. IV. In the appendix, 
we present details of the calculation. 

II. MINIMAL THEORETICAL MODEL 

To obtain the triple nonlinear momentum flux, (v r hu\\), we need to calculate the coherent 
part of fluctuations for the beat mode, and then use the two-scale direct interaction ap¬ 
proximation (TSDIA).— ^ In this way, the nonlinear parallel momentum flux can be written 
as 



( 3 ) 
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Here, the subscript (c) means the coherent component of the beat mode. v r = —ik y —( 

B 


is 


n 


the radial fluctuating E x B drift velocity. For simplicity, the adiabatic approximation— = 

n 0 


e0 . n 

— is used, so we have 


(c) 


n 0 


~t7 


Thus, we take mode widths large enough so that 


k^A 2 Vt he /(u}kV e i) > 1, with E the radial derivative of parallel wave number, A the mode 
width, Vthe the electron thermal velocity, and v e i the electron-ion collision frequency. k\\ = 0 
driven mode effects are neglected here. Now, the coherent parts of 0 and u\\ are required. In 
this work, we adopt three dimensional Hasegawa-Mima (H-M) equation with parallel flow 
compression which can be written as 

3 

— (p 2 s Vlcj) - 0) + p s u ci z X V0 • VVi0 - iu ;* n 0 = c,V||U||, (4) 


and the parallel momentum equation for cold ions 

3 \ 3 3 

— + u ciP 2 s z x V0 • V ) U|| - p s —{U\\}—ci> = -c fl V||0. (5) 

Here, we have used the standard normalization for electric potential fluctuation 0 = e0/T e , 
parallel velocity fluctuation u\\ = {tn/c a , with oj ci = eB/{rriiC ) the ion gyrofrequency, c s the 
ion acoustic velocity, and p s = the ion Larmor radius at the electron temperature. For 
the spatial scale, we consider two-scale approach, i.e., Vj_ = ik_L + d/dr, where kj_ denotes 
wave number of the fast spatial fluctuations, and d/dr describes modulation of the wave 
envelope, which occurs on a slowly varying spatial scale. oj m = k y p s c s /L n is the electron 
diamagnetic drift frequency with L n = — (d In n/dr) 1 density gradient scale length, and 
(U ||) is the mean parallel flow velocity. The last term on the RHS of Eq. (j3j) comes from ion 
parallel compression. In Eq. (J5]), the assumptions of isothermal electrons and u>k S> k\\{U\\) 
are used, and ion pressure gradient, Vy Pi, is absent due to the cold ion approximation. 

Taking the Fourier transformations of Eqs. fjTl)- (l5]) yields 

l 4+! rr|^^ + TTtrf“‘ = £ (6a) 

±rs ±rs k=k'+k" 

d d 

—'U/c — ikyp s -^(U\\)(f>k + ik\\c s (j)k — ^ ^ 

k=k'-\-k" 


(6b) 
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where the nonlinear terms are 

^ ri 


M *’ k '’ k " 2(1 + k\pl) 

d 


p 4 s {zx k'_L ■ k" (fcf - k'l)(t) kl (t) k n 


[ky{kf - k'l ) - x k' ± ■ k"] 


<9r 
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[KW-t - k ±) - 2k xZ x k ± ■ k±] 


(7a) 


M l,k',k" X k ± ' k l Pi {<t>k'U k » ~ U k '4>k") 


CO, 


T r\ ik v Ps ( V’k" ^ k ! r% 0fc 


2/Us 


d 


d 


dr 


dr 


to, 


a ik"p 2 s ( Uk"jr-4>k' — 4>k"-^-Uk' ) • (7b) 


d 


d 


dr 


dr 


d 2 


Here, the higher order terms related to slow spatial variation —have been neglected. 

dr z 

Eqs. (IbaH and (j6b|) can be expressed compactly in the form of a matrix as follows. 


d 


-g?&+H?4= y Kw 


( 8 ) 


with 


Vk = 


k=k'+k" 

4>k 

V*k 






r -] 

H = 

* i + /ci/| 

1 + 

= 

a b 


_f/c||C, -ikyPs — (f/||) 

0 


c 0 


and 


E M ‘. 


k'k" — 


k=k’+k" 


=fc'+fc" 

Mk,k', 


k=k'-\-k" 

The linear theory of this three dimensional Hasegawa-Mima system is clear. The disper¬ 
sion equation is as follows: 


. . iw.» . k V* 

Al = m ’“ ~ uiv + l TT' 

k?,c 2 

X 2 = ito k2 « -i- 11 —. 


(9a) 

(9b) 


The nonlinear terms are crucial to produce the coherent parts of the beat mode. By using 
the Markovian approximation, the nonlinear coupling terms can be written as 

E = - 7 lint + n,« + 2 


( 10 ) 
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where is the eddy-damping rate, and F kj0l is fast fluctuating force which does not con¬ 
tribute to coherent parts of the beat mode. Note that the nonlinear damping rate is larger 
than the frequency mismatch for strong turbulence. By diagonalization of the matrix H in 
Eq. (J8]), the coherent component of beat mode can be obtained as follows: 

% <0, (f)= [‘ «')<«». (ii) 

J — OO 

where the response function R^(t,t') is 

Rf(t, t) = rf S exp [( iu kS + 7 ^) (t' - t)] , (12) 


with 


a/31 

' k ~ 


a/3 2 
r k 


k\\C s 

1 

UJ * n 

k\\c s (l N k\p 2 s ) k\c 2 s (l + k\{%) 

^*n 

' k\c 2 s { 1 + k\p 2 s ) -fc,|C s ~ 
i ^*n ^*n 

—k\\c s (l + k\p 2 ) 




The details of calculation are presented in the Appendix A. 

Inserting the coherent component, Eq. (TTIT) into the nonlinear momentum flux, Eq. ([3]), 
then we need to calculate the forth order moment terms. By using the approximation 
of quasi-Gaussian statistics, the forth order moment can be decoupled into a product of 
quadratic moments, i.e., 


(Vkft)ril\t')g k ,,{f)ril,,{t')) = (pk'{t)p* k ,{t')){p kll (t)r]* k ,,[t')). (13) 


Here, within the Markovian approximation, the quadratic moments, or in another word, the 
two-time correlation function can be expressed by one-time correlation functions as 

(Vk*(t')v k (t)) = exp[iuj k (t’ -t)- 7j^| t' - t\]{ri%*(t)r)P k (t)). (14) 

Now, we have all the essentials for evaluation of the nonlinear momentum flux. 


III. NONLINEAR RESIDUAL STRESS AND COMPARISON 

In this section, we present the results for nonlinear momentum flux without showing the 
tedious calculations. The details of calculation can be found in Appendix B. Here, we write 
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the nonlinear parallel momentum flux again as follows: 

Ipif — t + (Vr^U^) + ■ 


(15) 


(c) (c) 

We need to substitute the coherent components (j) k into the first two terms and u k into 
the last term on RHS of Eq. to calculate the nonlinear momentum flux. The results of 
the first two nonlinear momentum flux terms can be written as 

n^ifi = (v^hu\\) + (v r n^u\\) 

= n 0 c 2 s $ {ik y - ik y)Ps[ 

7_, ;,// . J —OO 


k=k'+k" L 
_ „ NL d ( U II) 

- -"oXi -fir 


n 0 n 


NL,res 
r,||l 5 


(16) 


with the leading order nonlinear diffusivity is 


= 4^ ^ 1 + k\ p- 

k=k'+k" 


E 


r c 1 ^ci T T Ps a 

1 k’J-k" -j—gk"A-k’,k ", 
L J 


and the leading order nonlinear residual stress is 


y-[NL,res 1 2 

n r,||l = ~J n ° C s 


^ ^ Ik/Ikff 


^cl^ci 


L(i + fcip?) 


, A" 2 p, „L„_ lja _ s _ _ A" 2 

9kll-^kl ,kn T T T ' ^ t ^~c2^ci^y Ps9kl9kll r o 
L S L] Lj L s y L A 


k=k’+k" 

Here, r cl and r c 2 are triad interaction time for vorticity equation and parallel momentum 

equation, respectively. They can be estimated by the inverse of corresponding nonlinear 

damping rates, because the nonlinear damping rate is much larger than the frequency mis- 

1 1 Q 

match in strong turbulence. 7*. = 1 (j)A 2 is the fluctuation intensity, — = ——/& is the inten- 

L/ I k dr 

sity gradient scale length, L s is the magnetic shear scale length, A is the mode width, and 
other dimensionless parameters are Ak\k" — K ( K~ k ± + 2k 2)pi 9k = , 2 K ° sU 2 k v 

Wfc + lk ) NL 2 ) 

Note that Lj is positive for increasing intensity from inside to outside, and L s is positive 
for normal magnetic shear. Ak\k" > 0 for k’ ± ~ 7", and 9k is always positive. Therefore, 
the leading order nonlinear diffusivity satisfies xi L > 0 f° r increasing intensity profile in the 
edge regime. The sign of depends on the sign of L S1 and so is negative for normal 

magnetic shear. 












The other nonlinear momentum flux term can be written as 


TlNL 
LL r , ||2 


/--(c)v 

= {v r nu\ ) 


n 0 c 2 s 3? Y -itfyPs dt\R 2 ^%t,t')M^* k ^ k ,,(t')(p k/ (t)n k n(t)) 


= -noX 2 


k=k'-\-k" 
NL®(U ||) 


dr 


, yjNL.res 

+ n 0 n r ||2 - , 


( 17 ) 


where the leading order nonlinear diffusivity is 

X 2 -^PsC-s ^ ^ ^"c2 ^dJ-k'P-k" (,9k 1 9k"') Py P s ~j^~ 1 

k=k'+k" I 

and the leading order nonlinear residual stress is 

n JVL,res 1 2 T T 1 '2 2 Ps 

n r ,||2 2 ( ' s / ; ^dl^cJ-k'Ik"9k"k y P s ■ 

k=k'+k" s 

Here, the sign of \ 2 L is n °f dear. We can rewrite its expression in terms of symmetric k' 
and k", oc (g k / — g k ») ( k® — ky 2 ) , which is positive. Then, one can find that the sign of 
X 2 L is the same as that of x ± L , he., it is positive for increasing intensity profile. However, 
the sign of nonlinear residual stress, n^’ res is opposite to that of W Y p re ' s , i.e., Ii^p res is 
positive for normal magnetic shear. 

To campare with the usual Reynolds stress, we also calculate it quasilinearlly. 


= n 0 {v r u\\) 


3? Y n 0 c 2 (ikyp s (p* k u k ) 
j-j -QL,res 


k=k'-\-k" 

'll) 


= n 0 -x 


dr 


( 18 ) 


with quasilinear diffusivity 


and quasilinear residual stress 


X QL = p s c s Y hkIk > 


k=k’+k” 


A 2 


n °if’ r “ = -c? V h k — — I k . 

r,|| S Z-/ L v Lr 


k=k'-\-k" 


s^I 


Here, h k = 


klp s C s \Xk,NL2\ 

{ U k + Xk,NL2 ) 


Parallel symmetry breaking induced by fluctuation intensity 


gradient™ is used for the quasilinear residual stress. The quasilinear diffusivity is positive 
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TABLE I: Comparison of nonlinear and quasilinear results. 



diffusivity( / o s c s ) 

residual stress (c 2 ) 

(®r«||> 

e 1 / 2 

_ e 5/2 

(Vr^huu) + (v r h^Uu) 

e 3/2 

_ e 7/2 

{v r nu || ) 

e 5/2 

e 5/2 


definite. Different from the nonlinear residual stress, the sign of quasilinear residual stress 
depends on both Lj and L s . Yl®^' res is negative for increasing intensity profile and normal 
magnetic shear. 

Before comparing the nonlinear results with the quasilinear ones, we clarify the orderings 
of typical parameters that we will take in the following. The relative fluctuation amplitude 
from mixing length estimate, i.e., kg- ~ VJT use< F For the spatial scales, k x ~ k y ~ 

k_ l ~ 1/A, L n ~ L/ ~ L, and 


A 


~ ~ e are used, with e « 1 a small 

L J L s 

ordering parameter, which are consistent with TORPEX parameters.— For the temporal 


scale, normalized real frequency is order of 


. - T - e. The triad interaction time can 

^yCs ) 

be estimated as the inverse of nonlinear damping rate, because the frequency mismatch is 

much smaller than the nonlinear damping rate, as mentioned before. It was shown that the 

order of magnitude of the nonlinear damping rate for vorticity equation could be estimated 

as t'm ~ i , k^plujk-— Comparing the nonlinear terms in vorticity equations 

1 + k ± p s 

and the parallel momentum equation, one can divide 7 /"/ by a factor of k\p 2 to estimate 

lk 2 1 i- e -; lk 2 ~ Uk- The mean flow is one order smaller than the ion acoustic velocity, i.e., 

(U\\)/ C s e. Then we can estimate the order of magnitude of the nonlinear and quasilinear 

diffusivity and residual stress based on above orderings. The results are listed in Table 0] 

From Table Q] we can see that the nonlinear diffusivity is smaller than the quasilinear 

one. However, the leading order nonlinear residual stress is in the same order as quasi-linear 

residual stress, but with opposite sign for increasing fluctuation intensity profile. Note that 

the dominant contribution to the nonlinear residual stress comes from (v r fiu^) which is 

due to the coherent component of u\\. This is because the order of the nonlinear interaction 

2 

coefficient for u\\ is higher than that for 0*. by an order k\p 2 s ~ ~ e. Moreover, the term 

proportional to j^Ukr or 4^Uki, in M 2 k , k „ results in the leading order nonlinear residual stress. 
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The details of calculation can be found in Appendix B. The radial derivative of u k contains 
the radial derivative of k \\, due to the radial position dependence of k \\. We also note that an 
asymmetric parallel fluctuation spectrum is not necessary for non-zero, nonlinear residual 
stress due to the radial derivative of fcy. This is in contrast to the quasilinear residual stress, 
in which an asymmetric parallel fluctuation spectrum is required. For instance, the factor 

A 2 

j—j— in the quasilinear residual stress is induced by fluctuation intensity gradient symmetry 
breaking, which is order of e 2 . This is a partial reason why nonlinear residual stress can be 
comparable to the quasilinear one. Another reason is the relative fluctuation amplitude is 
large in strong turbulence, i.e., \4>k\ 2 ~ jt ~ e. This is different from weak turbulence, for 
which \(j> k \ 2 ~ e 2 . 


IV. SUMMARY AND DISCUSSION 

In the present work, we have derived the nonlinear parallel momentum flux for a three 
dimensional coupled drift waves and ion acoustic waves system. A Markovian approximation 
has been used for closure modelling. We estimate the triad interaction times for strong 
turbulence conditions. We also compared the nonlinear results with the quasilinear ones. 
It is shown that the nonlinear diffusivity is smaller than the quasilinear one. The leading 
order nonlinear residual stress can be comparable to quasilinear one based on the orderings 
we choose. This indicates that taking into account nonlinear wave-wave coupling effects 
on parallel intrinsic rotation is important. It is known that non-zero quasilinear residual 
stress requires symmetry breaking, such as fluctuation intensity gradient^ which we adopt 
in this work. However, in contrast to quasilinear theory, we find that an asymmetric parallel 
fluctuation spectrum is not required for a non-zero nonlinear residual stress. 

According to our theoretical results, the nonlinear residual stress is not negligible for 
strong turbulence. It supports the strong nonlinear parallel momentum flux induced by blob 
ejection which was observed in TORPEX expcrmiment.— In general, nonlinear momentum 
flux is also of potential relevance to tokamak edge region where similar blobs in L-rnode are 
also observed. Therefore, it may be needed to include the effects of nonlinear residual stress 
induced by blob ejection on intrinsic rotation in tokamak experiments. Recent experiments 
on ASDEX-U found the triple fluctuation term, v r hv po i is dominant as compared to the 
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poloidal Reynolds stress in turbulent poloidal momentum flux induced by ELM bursts during 
an H-rnode discharge.— In our future work, we will extend this work to focus on theoretical 
models of the nonlinear turbulent poloidal momentum transport. We also plan to compare 
nonlinear residual stress models with J-TEXT measurements of intrinsic torque and blob 
populations. 

Finally, we note that both resonant particle momentum diffusivity and residual stress 
can be calculated systematically to 0((pj.) in perturbation theory, for weak turbulence. The 
analysis follows Manheimer and Dupree^ and is similar to that for anomalous heating.— 
Though nominally higher order in perturbation theory, the results are not negligible, on 
account of differing resonant particle populations at different phase speeds. 
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Appendix A: Linear response function of three dimensional Hasegawa-Mima system 


The linearization of three dimensional Hasegawa-Mima system can be proceeded as fol¬ 


lowing. 



As 


r ~l 

H = 

\ + k\p 2 s 1 + k 2 ± p 2 s 

= 

a b 


ik\\c s 0 


c 0 


where, we have neglected ik y p s j^.(U\\). 


a — A b 
c —A 


0 =>- A 2 — aX — be = 0 . 


A 



2 
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Ai = i 


ca* 


1 + k\p 2 s U), 


1.2 2 u2 2 

A/ll Ay 11 

+ iJLA^x 2 = -iJLL. 


L U* 


(Al) 


The transformation matrixes are 


P = 


b b 


A 2 — a —b 

Ai — a A 2 — a 

’ 6(A 2 -A0 

a — Ai b 


By direct diagonalization of matrix H in Eq. flS]) and with the help of Eq. (ITU]) , we can 
get the equation, 


d 


dt 

Then, the above equation can be solved ss 


jgr/T + (A t ,„ + 7 “) C = 2+ F' ka . 


/q(c) 

Vk = 


dt'exp [(A fc , a + 7 k.a) (*' ~ *)] + F k,a) , 


(A2) 


(A3) 


where rf£ = F^ a = P a }F k 

Because F k) a is fast fluctuating force, which does not contribute to coherent parts of the 
beat mode, we can neglect it here. Multiplication of Eq. (lA3j) by matrix P on the left, we 
can get the coherent component of the beat mode 


“(<0 _ p '£( c ) 

Vk ~ F a ,pV k ■ 


(A4) 


The explicit expressions can be written as 

I dt'exp [{iuJki + Iki) (*' - t)] Pu [Pi\Ml k , k » (t 1 ) + P{ 2 M lk',At')\ 

J — oo 

+ [ dt'exp [(iu k2 + Vk 2 L ) (t' - *)] p i 2 [Pn M k,k>,k»tf) + p 22 lM k,k',k"( t ')] , 


$ 


U k ] = I- oo dt ’ eX P [i iUJ kl + lkl L ) ( t ' - t)] P 21 [Pu M k, k ,,k,,(t') + P 12 lM k,k',k"( t ')\ 

+ /-oo dt ' eX P [( iUJ k2 + 7 k2 L ) (*' ~ *)] P 22 [ P 2lM% tk , tk ,'(t) + M l,k’■ 

Then, we can rewrite them compactly in the form of matrix as Eqs. (HIT) and (1I2|) easily. The 
matrix in response functions are r^ 131 = P a \P ^, and r k p = P a 2 P 2 ( 3 1 respectively. 
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Appendix B: Calculation of nonlinear momentum flux 


More detailed calculation of nonlinear contribution is given this section. The linearization 
of Eq. (| 6 b|) can be written as 


^k 


k\\c s k'y p 


gw 

dr 


Uk + i r Jk,NL 2 


~4>k- 


(Bl) 


Because 042 resulted from correction of the parallel compression is much smaller than u k i, 
we only take uj k \ here, and omit the subscript 1 for simplicity in the following calculation. 
First, we present the result of Ilgp without tedious calculations. 


ngp = (vl c) nu\\) + (fynWy) 

= n 0 c 2 JR Y { ik y- ik 'y)ps[ 


k=k'+k" L 


+ 


Y ik v p * [ dt ' ex P[{- iuj ki + lki L ) -t)] ( M k%\k''( t ') ( Pk'(t)u k ''(t)) 

k=k'+k" d ~°° 

n 0 c 2 JR Y lk yPs [ dt'exp[(-iu kl + -f kl L ) (*' ~ t)] (—MW"(WW M W)) 
1 ui I ;.// J -00 ^ *n 


k=k'+k" 


1 j j & Y\\) 

^TlOpsCs / J 

k=k'+k" 


^~cl^ci Ps a 

9kf/~f~^-kf,k/f 


dr \ (1 + k]_p 2 s ) Lj 


~^~c2^ci 


a2 p k „ A " 2 A ' 2 + A " 2 

Ps ' ' ' 


hk/hk/r 0 fcy PsQkrQkrr 




'y ] Ik/Iku 

k=k'+k" 

A ' 2 A " 2 L, 


L s Lj L s Lj 

WlWci 


2 y rs - 


„ A " 2 p. ,V.. ,ja. 2 _ . A ” 2 

(■\ j2 9 \9kr/Siki : kri , , T + A r c 2UJ c iPy PsQktQklt 

ft + ^±Ps) J- 1 ! J-'s 


L 2 


^~c2^ci~ 


NL ~ 11/ 1 TTivi^.res 

= + n ° u r,\\{ > 


L 2 L 2 


h k ,h kM + p s p fc „ - 3fc„ p s g k / 9 k// 


d{U\\) 


NL.res 


dr 


(B2) 


with 


xf L = 


^ ^ Ik/Iku % 

-U\un v 


/c=fc , +/c // 

L 


WlWcz Ps 

(1 + /cip 2 ) L/ 


Qk/lA-kr, 


kff 


7~c2^ci~ 


A " 2 A ' 2 + A " 2 

‘IPy PsPkH p k 


s^I 


l q l 


S^I 


hk,h k „ g PsPkrPkrr 
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and 


yrNL,res _ 1 2 

nil 1 ~ _ 4 C * 


^ ' Ik/Iktl 
k=k'+k" 

A' 2 A" 2 L„ 


L(1 + k\pl) 


A" 2 p s L n 2 2 A" 2 

9kn^kt,kff~^ j j H - ^~c2^ci^y Ps9kf9kff~ 


L s Lj Lj 


LI 


~T c 2C0 ci~ 


L 2 ,L ? L, 


hk/foktr ky p s Pk/f Ps9kf9kff 


y 2 p 2 s± 


Here, the dimensionless parameters are A k ^ k n = /c" 2 (A/J 2 — A;" 2 + 2k' 2 ) p A s , g k = 


k 2 p s c s ^ k:NL2 k 2 c 2 

h k = —o-v and Pk = 


( W fc + lk,NL2) 


order. 


(^fc + 7fc,ATL2) 


+ 7fc,JVL 2 ) 


. The first lines in xf L and n A ,^’ res are leading 


The calculation of n^£ is written as follows. 


n^|f 2 = (tv-rutf^) 

= n 0 c^3fj ^ -ikyPs f dt\Rl! 3 *(t } t')Ml!* k , tk ,,(t')(j) k/ (t)nk-(t)) 


k=k'-\-k" 


= n 0 c 2 Jk Y ikyPs j dt , exp[(-iu; k2 + 'r k2 L ) (t'-t)] (M k * k , tk/I (t')(t)k'(t)n k "(t)) 

k=k’+k" J _0 ° 

+ n 0 c 2 dt Y ~ ik yP s [ dt’exp [(-iuj kl + j kl L ) (t' - t)] ( 

k=k'+k" d-oo 

+ n 0 c 2 s ^t Y -ikyPs [ dt'exp [{-iu k2 + 7k> L ) _ 0 ] ( 

7 _ 1 un J —OO 


k\\c s (l + k 2 ± p 2 ) u 

k,k',k"V / )$k'\t)n k "(t)) 


„2„2(-t 1 7,2 „2t 




k=k'+k" 


-\no<?s& Y1 T ^dkyp 2 s (<l>l,(t)<i> k ,(t))(p a ^^ n k „{t)) - (ul,(t)<j> kl (t))(p s d ^^ n k ,,(t)) 


k=k'-\-k" 


dr 


dr 


+ 1 -n 0 clX V TclWd hAi±l±dh« (fc« _ *42 + 2k?) pM,(t)Mt))(p,^f^-n u (t)) 

k=k'+k " CU * r ’ r 


2 

1 


— -n 0 c 
2 ° ' 


® V 

k=k , -\-k" 


^~c2^ci~ 


^*72 

fc 2 c g(l + k\ p 2 s ) 

ujL 


dii* Hi 

K 2 pI W)Mt))(Ps-^^n kn {t)) 


dr 


- (uUt)MmP,‘LL^-n u (.t)) 

As we can see, we have written Ip£ into three components, and we will calculate each of 


(B3) 
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them one by one as following. 

du* k „(t) 




k=k'+k" 


dr 


n k n{t)) 


E 

k=k'-\-k" 

E 

k=k’+k" 


Ps 


, x" p & 


a < Vl) K». 


Co dr 


Ukn - *7^/2 


Ps 


dcf) 


dr 


k "<t>k„ ) ~ 


Ps J kyPs 


■ NL (fiknfikir) 

Vk» - *7fc/,2 


1 Kpsl 


kir 


Ps 


9(U\\) 


2 w fc// - 


2 — + 

L s Li dr 


(B4) 


X 

Here, k\\ = k y — is used. L s is the magnetic shear scale length, x = r 0 — r, where r 0 is 

L o 


the radial location of resonant surface. So, the second term in the second line comes from 
d k 

— k\\ = -Note that parallel asymmetric fluctuation spectrum is not required for the 

dr L s 

second term. 


E _ e±LEL )Kpi 
E ( U u(t)<t>u{t)) = ( ~ L ‘ .L-fl.NL - 


fc=/c'+fc" 


k=k’+k" 


Ps 


Uk, - *7** 


E 

fc=/c , +fc ,/ 


P>s4/ (c A' 2 a<[/||> 


Ukt-ilhn \PsL s L / 


+ 


<9r 


(B5) 


According to Eqs 
1 


and (IB5I) we can get the hrst component of nE4 i.e., Eq. ()B3j) . 


~2 n ° Ci 


3? ^ ^ r c2 C0ciky p 


r, ||2’ 
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dr 
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Lo Lr dr 
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+ 


dr 


2L, 


= c s 


^ ^ 7"c2^ci4/4// 


k=k’+k” 


,/2 


Ps / \ Ps //22^ (^11) , Ps j/2 A 

■ 7 “ \9ki ~ 9kn ) —P s —w-K - 7 —Pfc, 

iv/ Cs Or Li I X/gXiJ 
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Now, we turn to calculation of the second component of IlN^. 


\no<?s& ^ k]Ml + k±P ^ k y { k ± - k ± 2 + 2 k 'x) Pa(^fc/(*)^(*))(P« ^^ t ^ fc//(*)) 

k=k'+k" U * n 


dr 


o x Ln ~t~ x") (1 -\- k~i p 2 ) , , . / w i dcbT/ili) , .. 

n 0 c^ r c ia; ci -—- -A kl ^ kl ((j) kl (t)^ k ,(t))(p s — ■^~ n kn{ t )) 


k=k'+k" 


= 


E 


Lsps 

_ _ L n (l + A^p 2 ) A ' 2 

'cl ^ ci -™-k//,kr ^kt-*-kih 

s ^ j 


dr 


Noting that we have used the relation 3? (x^p, 

above. 


k=k'-\-k" 


dr 


n kll (t )) = 0 in the calculation 


By the same token, we can get the last component of IlN£ as follows: 


"2 " 0C; 


;* V 


T&Wd- 


Bcl{l + k 2 ± fj 


ui 


k ,2 n 2 

^y Ps 
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,2 

SPS 
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<9?/* 1/1 
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\n 0 c 2 s ^ T c2 u ci I kl I kll j±( 1 + P 2 Pa) j ^ 2 (A /2 + A" 2 ) 
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Kps 


Cs ^ Ps 


dW' 


u k n - ?7fc//2 v L * 2L I d' r 


+ 


*>. /c, A“ a(f/||>\ p. 


+ 


Vki-ilka \psLsLi ’ <9r y 2L 7 

1 2 T T L n( l + ^±P 2 ) 

_ o n 0 C a 2^ T c2 UJ ct I k ,I k jr 


- k 2 A 


/2 A /2 y 


Kps A" 2 3c, 


L 2 L 


k=k'+k" 

A p ‘ A “ 


L 2 


k' 2 (A' 2 + A" 2 ) 


^kn l lkn2 

p s P, B(U t |> 
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c, 2L/ <9r 


‘2L, Lj 


2 Pfc r 

2 E, 


a" 3ps k’ 2 A' 2 A " 2 
9k T K y ^ r 2 f ’ 


(B6) 


Here we have used the relation: 


so /„»„ du* k „(t) ^ _ „ T A" 2 3p s 

J? / ^ Ps TOjb,/(t)> Qfclkh 


k=k'+k" 


dr 


L s 2Lj 
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Finally, we combine all the components of II^p to obtain 


■^r,||2 , TloPs^s ^ ^ d^k/Ikir (Pkr Qkll ) 


k=k'+k" 


by] + y 1 2 3 4 5 6 ( a' 2 + a" 2 )!” (1 + ilft) 


H - ^ ^ ^~c2^cilk/lkff ^ Qkr/ky Ps 


k=k r +k" 
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||) 


2L 2 


(1 + 


LI 

Ls{ 1 + fcjpg) 

L n k 2 ± p 2 s 


Akir,kt T 9 

L 2 


Ps d(U\\) 
Lj dr 

A' 2 


A' 2 


Ps 


L2 (1 + fc 2 p 2 ) L 2 L, 


= -^0*2 


(9r 


■ -i-rTVL.T’es 

+ n 0 n r||2 , 


(B7) 


with 


X2 L = \psC s 


^ T::2 : ~^ci^-kl kklf (Qkl Qkll) 


k=k'+k" 


k ,2 p\ + k ,2 ( A' 2 + A" 2 


L 2 


Ps_ 

Lj ’ 


and 


yjNL,res 2 

n r,p = 2 C * 


^ ^ ^~c2^ci^k/^kn A Qkf/k, 


'2 2 Ps_ 
y Vs t 


k=k'-\-k" 

-2g kr ,k' 2 (A' 2 + A" 2 ) £ + <fc,A*pj 


2L2 


(l + k 2 ± p 2 s ) 


Lnk“^_p 2 k,, ' k ’ Lj 


A ' 2 


Ps 


LI (1 + kip]) Lj L s \ 


Note that the first terms of xl L and IF V |^’ res are leading order. We have represented r c i 


r, \\2 

2 „2 


kip. 


■T c i which comes from the comparison between 


by r C 2 by using the relation r c 2 = , , 0 „x 

(1 + klpj) 

nonlinear terms in vorticity equations and those in parallel momentum equation. 
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